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Abstract. For every nonholonomic manifold, i.e., manifold with nonintegrable distribu- 
tion the analog of the Ricmann tensor is introduced. It is calculated here for the contact 
and Engel structures: for the contact structure it vanishes (another proof of Darboux's 
canonical form); for the Engel distribution the target space of the tensor is of dimension 2. 
In particular, the Lie algebra preserving the Engel distribution is described. 

The tensors introduced are interpreted as modifications of the Spencer cohomology and, as 
such, provide with a new way to solve partial differential equations. Goldschmidt's criteria 
for formal integrability (vanishing of certain Spencer cohomology) are only applicable to 
"one half" of all differential equations, the ones whose symmetries are induced by point 
transformations. Lie's theorem says that the "other half" consists of differential equations 
whose symmetries are induced by contact transformations. Therefore, we can now extend 
Goldschmidt's criteria for formal integrability to all differential equations. 



Introduction 

Let Wff be the space of germs of irreducible k- dimensional distributions at the origin of 
". The orbits of DiffIR n -action on W£ are usually of infinite codimension. As shown in 



VG1|| , there are 3 types of exceptions: (1) Wff is a point, (2) codimension 1 distributions 
determined as zero set of a = for the forms a of maximal class (recall that the class of the 
form is the invariant equal to the least number of variables actually involved) 



a 



dt + (.Pidqi — qidpi) for n = 2r + 1 

l<j<r 

dt+ Yl {Pidqi - qidpi) for n = 2r 

KKr-l 



called (usually, for n odd) a contact structure (see |Z|) and (3) the Engel distribution locally 
given by the system of Pfaff equations a\ = 0, a 2 = 0, where 

q;i = dx4 — x%dxi, c*2 = dx3 — x<idx\. (E) 

Any 4-dimensional manifold with the Engel distribution (E) will be called an Engel manifold. 

Both contact and Engel manifolds are examples of nonholonomic manifolds (i.e., mani- 
folds with nonintegrable distributions). Recall that the term nonholonomic was coined by 



H. Hertz [H] in his attempts to geometrically describe motions in such a way as to eliminate 



the notion of "force". For a historical review see |[VG2j| ; see also a very interesting paper 



by Vershik [0] with first rigorous mathematical formulations of nonholonomic geometry and 
indications to various possible applications in various, sometimes unexpected, areas (like 
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optimal control or macro-economics). In Vershik summarizes about 100 years of studies 
of nonholonomic geometry (Hertz, Caratheodory, Vranceanu, Wagner, Schouten, Faddeev, 
Griffiths, Godbillon, to name a few; to his list I'd add that MathSciNet returns about 700 
entries for nonholonomic, and for synonyms: 70 more for sub-Riemannian,1700 for Finsler, 
several dozens for "cat's problem" and autoparallel; physics have their own simialr lists with 
thousands of entries because there seems to be more, actually, nonholonomic dynamical 
systems than holonomic ones; finally supergravity is also a nonholonomic structure, albeit 
on supermanifolds — ca 1800 more entries). In |V[] Vershik also declares that "it is proba- 
bly impossible to construct an analog of the Riemann tensor for the general nonholonomic 
manifold" though in some particular cases of small dimension they had been constructed 
(Martinet) (and, one can add now, in papers on supergravity). 

Here I introduce the invariants that characterize "nonflatness" of nonholonomic manifolds, 
i.e., positively solve this old problem. This definition appeared around 1989 but publication 
of this paper was delayed and the definition first (as far as I know) appeared in contetext 
of supergravity in ||LP|| . For application to SUGRA (supergravity equations) and other 
instances where the N-curvature tensors (the left hand side of SUGRA) are calculated see 
I GLT1 , jLPSl , jLSVfl. 

The general definition is illustrated with contact and Engel manifolds. As expected, for 
the contact structure these invariants vanish, thus giving one more proof of the Darboux 
theorem on the canonical form of the contact form; one can show that they vanish for n even 
as well. The Engel manifolds are more rigid: there could be nonflat ones among them. 

Observe that our invariants, though natural analogs of the curvature and torsion tensors, 
do not coincide on nonholonomic manifolds with the classical ones and bearing the same 
name. Indeed, by definition on any nonholonomic manifold, there is a nonzero classical 
torsion while every contact manifold is flat in our sense. To avoid confusion, we should 
always add adjective "nonholonomic" for the invariants introduced below. Since this is too 
long, we will briefly say N-torsion, N-curvature, etc. 

For manifolds with a G-structure, the prototypes of the invariants we consider are well- 
known: these are structure functions of the G-structure |St[ ; an important particular case is 
G = 0(n) when the structure function is called the Riemann tensor. 

A theorem of Lie states that there are two types of differential equations, considered as 
submanifolds in the jet bundles: (1) those whose symmetries are induced by point trans- 
formations and (2) those whose symmetries are induced by contact transformations. For 
equations of the first type Goldschmidt gave criteria for formal integrability | [BCG| . The 
notions introduced below — generalizations of Spencer cohomology for nonholonomic case 
- provide with a straightforward generalization of Goldschmidt 's criteria for formal inte- 
grability to the second half of differential equations. 

Finally, observe that "non-flatness" of distributions with non-linear constraints (such as 
fields of cones fA|, foreseen by Vershik in 0, or fields of spheres, like the ones constructed by 
any driver who switched the cruise control ON, can also be tested by invariants I introduce; 
in the case of the fields of solids, however, the Lie algebras involved are of infinite dimension 
(they correspond to "curved Grassmannians" of submanifolds, rather than linear subspaces) 
and feasibility to obtain the final result becomes very doubtful. 

The waves in plasma is another example of a nonholonomic system with infinite dimen- 
sional algebras involved. This is a pity, since the curvature tensor is known to be responsible 
for stability of the dynamical system expressed by geodesies; so, in principle, we could just 
have computed the N-curvature in this case and settle an intriguing question: "is cold fusion 
possible or not?". 
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§1. Structure functions 

Let M be a manifold of dimension n over a field K. Let F(M) be the frame bundle over 
M, i.e., the principal GL(n; K)-bundle. Let G C GL(n;K) be a Lie group. A G-structure 
on M is a reduction of the principal GL(n; K)-bundle to the principal G-bundle. Another 
formulation seems to be clearer: if one can select transition functions from one coordinate 
patch to another so that they belong to G for every intersection pair of patches, then we 
have a G-structure on M. 

The simplest G-structure is the flat G-structure defined as follows. Let V be K n with a 
fixed frame. Fix a frame in V. Having identified V with T V V via parallel translation along 
v, we get a fixed frame in every T V V. The flat G-structure is the bundle over V whose fiber 
over v E V consists of all frames obtained from the fixed one under the G-action. 

In textbooks on differential geometry (e.g., in [St]) the obstructions to identification of 



the k-th infinitesimal neighborhood of a point m E M on a manifold M with G-structure 
with that of a point of the manifold V with the above flat G-structure are called structure 
functions of order k. It is also explained there that such an identification is possible provided 
all structure functions of lesser orders vanish. 

To formulate the precise statement, we need some definitions. Let S l denote the operator 
of the i-th symmetric power. Set 

_! = T m M, 0o = = Lie(G); 

for any i > and k — 1, . . . , i set: 

Qi = {X E Hom(g_ 1 ,0 i _ 1 ) | X(v x )(v 2 , v 3 , v t ) = X{v 2 ){v x , v 3 , u 4 ) for any v k E g_i}. 

Now, set (fl_i,fl )* = © 8i- 

i>-l 

Suppose that the g - m °dule g_i is faithful. Then, clearly, 

(fl-ij 0o)* C t)cct(n) = Oct K[[xi, . . . , x n }}, where n = dim g_i. 

It is subject to an easy verification that the Lie algebra structure on oect(n) induces same 
on (fl-i,flo)*- The Lie algebra (fl-i,go)* wn l be called the Cartan's prolong (the result of 
Cartan's prolongation) of the pair (g_i,go). 

Let A 1 be the operator of the i-th. exterior power; set C,£_^ g( s = gfc-s®A s (gI 1 ); we usually 
drop the subscript of G^ s i 0o ^ or indicate only go- Define the differential d s : C k,s — > Qk,s+i 
setting for any v x , . . . ,v s+x E V (as usual, the slot with the hatted variable is ignored): 

(5 s f)(v 1: ... ,v s+1 ) = ^(-l) l [/(fi,... ,v a+1 _i,... ,v a+1 ),v a+1 -i]. 

i 

As expected, S s S s+ i = 0. The homology of this bicomplex is called Spencer cohomology of 
the pair (g_i,g ) and denoted by H^ s _ i So ^. 

Proposition . ([f5tj|) Structure functions constitute the space of the (k,2)-th Spencer coho- 
mology of the Cartan prolongation of the pair (T m M, g) where g = Lie(G). 

The Riemann tensor on an n-dimensional Riemannian manifold is the structure function 
from H? 2 , \ \ . 

(fl-i,o(n))» 



Remark . (cf. 



k 



^Sl 1 , fc ). = ^ 5 (fl-i;(fl-i.flo).). 



4 



DIMITRY LEITES 



This remark considerably simplifies calculations, in particular, if the Lie algebra (g-i, flo)* 
is simple and finite dimensional: we can apply the Borel-Weil-Bott- . . . (BWB) theorem. 
In the nonholonomic case considered in what follows the remark provides with a compact 
definition^ of structure functions. We can recover the bigrading (explicitly preserved in 
0) at any moment but to work with just one grading is much simpler, as with generating 
functions of series instead of individual terms. Moreover, there are a few general theorems 
that help to calculate Lie algebra cohomology and no theorems on Spencer cohomology. 



§2. Structure functions of nonholonomic manifolds 

To embrace contact-like structures we have to slightly generalize the notion of Cartan 
prolongation: the tangent space to a point of a nonholonomic manifold is naturally rigged 
with a nilpotent Lie algebra structure (e.g., for any odd dimensional contact manifolds this 
is the Heisenberg algebra), cf. HVG ]. 



2.1. Nonholonomic manifolds. ( ||VG1|| , |[VG2;1 ) Let M n be an n- dimensional manifold 
with a nonintegrable distribution D. Let 

D = D 1 cD 2 cD 3 -- - cD d 

be the sequence of strict inclusions, where the fiber of Di at a point x G M is 

(here [Di, -Dj_i] = {[X, Y] | X G Z?i, Y G -Dj-i}) an d d is the least number such that 

D d \ x \J [Di,D d ]\ x = D d \ x . 

In case D d = TM the distribution is called completely nonholonomic. The number d = d(M) 
is called the nonholonomicity degree. A manifold M with a distribution D on it will be 
referred to as nonholonomic one if d(M) ^ 1. 

Here we will consider nonholonomic manifolds only when they are completely nonholo- 
nomic. This does not diminish generality, actually, since if the distribution in question is 
not completely nonholonomic, we can take the integral manifold corresponding to it; the 
restriction of the same distribution onto this integral manifold is completely nonholonomic, 
all right. If, nevertheless, we have to, for some reason, consider the general case it is easy 
mutatis mutandis; it corresponds to nontransitive Shchepochkina prolongs, cf. below, and 
the general machinery is easily adjusted to them). 

Let Ui(x) = dim Di\ x . The distribution D is called a regular one if all the dimensions ni 
are constants on M. In what follows we will only consider regular distributions. This is also 
not a restriction if everything is smooth. 

The tangent bundle over a nonholonomic manifold (M, D) can be naturally endowed with 
a bundle of nilpotent Lie algebras as follows. Let L be the Lie algebra of vector fields 
preserving the distribution in a neighborhood of some point x G M. Naturally, L is filtered 
and let g be the associated graded Lie algebra. Set 

0- = © 0*. 

-d<i<-l 



Clearly, d, the depth of g_, coinsides with the nonholonomicity degree of D, see 

An analog of G-structure for nonholonomic manifolds is defined to be the pair (g_,go); 
where Qo is a subalgebra of the Lie subalgebra which preserves the Z-grading of g_, i.e., 



^f. with the problems encountered in the pioneer paper Q, where the case d — 2 for Lie algebras was 
independently introduced. 
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go C (9er 0-)o- If 0o is not indicated explicitely we assume that q = (Det 0-)o, i-e., is the 
largest one. 

Thus, we see that with every nonholonomic manifold (M, D) a natural G-structure is 
associated: Lie(G) = (Det 0_)o- But the structure functions of this G-structure do not 
reflect the nonholonomic nature of M. Indeed, when we compute the structure functions (be 
they "traditional" H,^., or introduced above H s (g_i, (fl-i,flo)*)) we use the fact that 
0_i is a commutative Lie algebra (partial derivatives commute). 

To take the nonholonomic nature of M into account we need an analog of the above 
Proposition for the case when the natural basis of the tangent space consists not of partial 
derivatives but rather of covariant derivatives corresponding to the connection determined 
by the same Pfaff forms whose vanishing determines the distribution. Therefore, instead of 
0_i we have a nilpotent algebra 0_. 

So we have to define 

(1) the simplest nonholonomic structure — the "flat" one, and 

(2) the analog of Cartan prolong when g_i is replaced with g__, 

(3) (we have already understood) what is the analog of O if it is not mentioned as is 
usually the case (only the distribution is given), but it remains to figure out 

(4) what is the analog of (0-i,0o)*), and, finally, 

(5) what is the analog of H?' z „ % . 

Let us answer these remaining questions. 

Fix a frame F in K n . Having identified T V V with V by means of the transformation that 
preserves D, we fa a frame — the image of F — in each T V V. The nonholonomic structure 
of the nonholonomic manifold (K n , D) is the pair of bundles: the frame bundle and the 
distribution D. The nonholonomic structure of (K n ,D) will be considered to be flat if the 
frames over v are obtained from the fixed one by means of the G-action, where G is a group 
whose Lie algebra is (Det 0-)o- 



2.2. Structure functions of a nonholonomic manifold. Given a pair (g_,0o) as above, 
define its k-th Shchepochkina prolong (first published by I. Shchepochkina for Lie su- 
peralgebras; see also more accessible ||Sh5|| ) for k > to be: 

g k = (i(S*(g^y <g> O ) n j(S*te-T ® 0-))ft, 

where the subscript singles out the component of degree k and where i : S ,fe+1 (g_i)*(S>0-i — ► 
-i)* ® 0!_i ® 0-i is a natural embedding and j : S k (g-i)* ® 0o — ► 5' fc (0-i)* ® 0* i <8> 0-i 
a natural map. 

Similarly to the case when 0_ is commutative, define (0_,0o)* to be © 0^; then, as is 

k>~d 

easy to verify, (0-,0o)* is a Lie algebra. By the same reasons as for the manifolds with G- 
structure (0), the space i? 2 (0_; (0_, 0o)*) is the space of structure functions — obstructions 
to the identification of the infinitesimal neighborhood of a point on the manifold with a 
nonholonomic structure (given by 0_ and 0o) with same on a flat nonholonomic manifold 
with the same 0_ and O . 

The space H 2 (g_; (0_,0o)*) naturally splits into homogeneous components whose degree 
will be called the order of the structure function; the minimal order of structure function 
is equal to 2 — d. As in the case of a commutative 0_ = 0_i, the structure functions of 
order k can be interpreted as obstructions to flatness of the nonholonomic manifold with 
the (0_, 0o)-structure provided the obstructions of lesser orders vanish. Observe that for a 
nonholonomic manifold the order of structure functions is no more in direct relation with 
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the numbers of the infinitesimal neighborhoods of the points we wish to identify: distinct 
partial derivatives bear different "weights" . 

When I told this to A. Beilinson in 1989 he winced and, as I got this far, he said that 
though he conceeded the idea is of interest, he doubted it can work in the general case: all 
my constructions are graded, whereas with every distribution only filtered algebras g_ are 
naturally associated. This incredulity is justified but fortunately, different nitrations on g_ 
with the same graded $j_ are governed precisely by the structure functions of lesser orders 
and all of them must vanish for the N-curvature to be well-defined. So, given a filtered 0_, 
just take the associated graded algebra and proceed as indicated above. This indicates a 
way to study the stable normal forms of Pfaff equations listed in B) . 



§3. The Lie algebra e of Engel vector fields 

Let us describe the Lie algebra e of vector fields that preserve the Engel distribution. Since 
each of the forms oti and a 2 is a familiar contact form, we know enough to demonstrate (we 
leave it as an easy exercise for the reader) the following statement. 



Lemma 

grees): 



A weight basis of the Z-graded Lie algebra e is as follows (indicated are the de- 



-3 


-2 


-1 





n > 1 


d, 


d 3 + xid 4 


dx 


xidi + x 2 d 2 + 2x 3 d 3 + 3x 4 <9 4 


x n+lQ + *i Q + *i Q 
J- z n+2 ° n+3 ^ 






d 2 + xtd 3 + 


£i<9i — x 2 d 2 + x 4 <9 4 
%id 2 + ^d 3 + ^d 4 



The structure of the Engel algebra e seem to have never been described; so let us do it. 
First, introduce a shorthand for some elements of e: 



D 3 
E 
H 



d 3 + £i<9 4 ; 

xidi + x 2 d 2 + 2x 3 d 3 + 3a; 4 <9 4 ; 
xidi - x 2 d 2 + x 4 9 4 ; 



x'i ' d 2 + - J T7rd 3 + - j tt9a for n > — 1. 

J- z n+2 ° n+3 * — 



Now, it is clear that e is solvable with 

e/[e,e] = SpaniE^H); e (1) /[e (1) , e (1) ] 



Span(di), where c^ 1 -* = [c, e]. 



plays for e the same role the Heisenberg algebra plays for the contact 

|G1 | where it also bears 



The negative part, e 
algebra. This negative part was, of course calculated earlier, see, e.g. 
the name of Engel algebra; perhaps one should say the "small" and the "total" ones. It seems 
very strange that the total algebra had never beed desribed; I am sure such a description is 
buried in works of classics. 

In review to [|G1|| [MR 96k:58010] it is written "There is a Darboux theorem for Engel 
structures: any two are locally diffeomorphic" . Vershik also quoted this statement as a 
well-known. But I can not place it. Indeed, by Darboux theorem for the contact structure 
I understand the statement that says that the contact form can be locally reduced to a 
canonical form , cf. |ZJ. (In contradistinction, a related statement that the nondegenerate 
differential 2- form uj can only be reduced to a canonical form if and only if dw — 0, cf. 0, 
so the statement on reducibility of the nondegenerate closed 2-form, sometimes encountered 
in textbooks on Differential Geometry, is a tautology.) 

Concerning Engel structures, there is, for example, paper \U2\ with examples of several 
non-equivalent Engel structures on M 4 . Doesn't it contradict the statement on an analog of 
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Darboux theorem? So let us calculate the N-curvature of the contact and Engel structures: 
to find out if these invariants vanish, hence, there are analogs of Darboux theorem for them. 

§4. Structure functions of the contact and Engel structures 

Eli Cartan worked before "(co) homology" was defined, still, he introduced a notion that 
helps to calculate some of these cohomologies. This notion is involutive Lie algebra of 
vector fields. Serre gave a nice reformulation of this notion (Pt|); thanks to Serre's state- 
ment for the case when (fl-,flo)* is a simple vectorial Lie algebra instead of computing 
H 2 (Q-i, (fl-i,0o)*) — ®-^(g 2 i 0o ).' ^ on ^ suffices to compute the summand with the least 

k = 1: only this term might be nonzero because (g_, g )* is involutive. This really helps, see 
examples in ||LPS|| . 



Regrettably, Serre's theorem is only formulated so far for algebras of depth 1. So I cal- 
culated H 2 (q-; 0o)*) by hands; later the result was confirmed by independent computer 
calculations up to convincingly high degree: 

4.1. Theorem . The space of nonholonomic structure functions on any Engel manifold is 
of dimension 2. All structure functions are of order 2, their representatives are: 

<9* A X*_ x <g> H + <9* A X*_ x ® E and d* 4 AD;® <9 4 . 

The first of these cycles is "pure" in the sence that it does not represent its cohomology class 
modulo any coboundary, the coboundaries vanish. The other one is mixed: one can add any 
coboundary to it. 



The following is another form of the classical theorem (see |LF|): 



4.2. Theorem. (Darboux theorem) The space of nonholonomic structure functions on any 
contact manifold vanishes. 
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